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matlcal  formulation  was  accomplished  by  Prof.  T.W.  Ting  of  the  University 
of  Illinois  while  working  as  a  visiting  scientist  at  the  Applied  Mathe- 
matics  Laboratory  of  the  Aerospace  Research  Laboratories,  Wright-Patterson 
Air  Force  Base,  Ohio.  A.  F.  Grandt  was  the  Project  Engineer  for  the  Air 
Force  Materials  Laboratory.  The  report  covers  research  conducted  under 
Project  No.  2307,  Task  No.  2307P1 ,  Work  Unit  No.  2307P102.  This  report 
covers  work  conducted  from  January  1975  through  December  1978. 
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SECTION  I 
INTRODUCTION 

Since  mechanically  fastened  joints  represent  one  of  the  most  common 
failure  sources  in  aircraft  structures  (Reference  1),  considerable  re¬ 
search  has  been  directed  toward  developing  fatigue  resistent  fastening 
systems.  One  common  method  to  extend  the  fatigue  life  is  to  introduce 
a  controlled  residual  stress  field  around  the  fastener  hole  by  means  of 
an  interference  fit  fastener  (Reference  2)  or  by  pre-expanding  the  hole 
with  an  oversize  mandrel  (Reference  3).  This  latter  process,  commonly 
known  as  coldworking,  differs  from  interference  fit  applications  mainly 
in  that  the  hole  is  allowed  to  radially  unload  prior  to  inserting  the 
fastener.  This  unloading  results  in  a  large  compressive  stress  field 
next  to  the  edge  of  the  hole  which  improves  fatigue  life.  Coldworking 
precracked  fastener  holes  has  been  shown,  for  example,  to  increase  the 
cyclic  life  by  two  orders  of  magnitude  in  some  cases  (Reference  4). 


In  order  to  analytically  predict  the  service  life  of  coldworked 
holes,  it  Is  necessary  to  examine  the  residual  stress  field  caused  by 
the  pre-expansion  process.  An  early  study  of  the  plastic  deformation  of 
holes  In  flat  rings  and  disks  is  discussed  in  Chapter  33  of  Reference  5. 
Although  the  applications  discussed  In  that  report  center  about  the 
practice  of  Interference  fitting  small  tubes  into  industrial  boilers  and 
condensors,  their  form  of  solutions  for  stresses  in  the  plastic  region 
are  employed  here  for  the  coldworked  fastener  hole  problem.  Several 
authors  have  examined  the  stress  fields  around  interference  fit  fasteners 
by  both  experimental  (References  2,  6-7)  and  analytical  or  numerical 
procedures  (References  8-10).  However,  analysis  of  unloading  which 
characterizes  the  coldworking  process  (References  11,  12)  has  received 
less  attention  In  the  published  literature.  The  problem  formulation 
and  mathematical  analysis,  which  provide  a  basis  for  the  results  pre¬ 
sented  here,  are  reported  in  References  13  and  16. 
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The  major  content  in  this  report  is  in  two  sections,  Elastic-Plastic 
Analysis  and  Discussion  of  Results.  In  the  Elastic-Plastic  Analysis  section, 
deformation  theory,  radial  symmetry  and  plane  stress  assumptions  are  used 
to  model  the  coldworking  process.  Comprehensive  results  for  finite  ring 
geometries  and  arbitrary  material  properties  are  presented  in  graphical 
form.  In  the  Discussion  of  Results  section,  the  features  of  the  graphs 
which  seem  most  significant  from  a  designer's  viewpoint  are  emphasized. 

The  graphs  Indicate  that  the  region  of  yielding  during  pre-expansion  char¬ 
acterizes  the  residual  hoop  stresses  which  presumably  inhibit  crack  growth. 

A  relation  between  this  region  of  yielding  and  the  mandrel  interference 
leads  to  an  analytical  optimum  interference.  A  comparison  between  analy¬ 
tical  optimums  and  experimentally  determined  mandrel  interferences  which 
produce  maximum  fatigue  life  shows  very  good  agreement.  Thus,  it  is 
suggested  that  the  cause  of  the  beneficial  effects  of  coldworking  is  cor¬ 
rectly  characterized  and  the  analytical  optimum  interference  corresponds 
to  the  mandrel  Interference  that  leads  to  maximum  fatigue  life.  The  major 
ideas  are  Included  in  the  summary. 
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SECTION  II 

ELASTIC-PLASTIC  ANALYSIS 

1.  MATHEMATICAL  MODEL 

The  geometry  for  which  the  coldworking  process  is  modeled  is  depicted 
schematically  in  Figure  1.  We  consider  a  hole  of  radius  a  in  a  circular 
sheet  of  elastic-perfectly  plastic  material  with  outer  radius  b.  The 
coldworking  process  involves  radially  expanding  the  hole  with  a  mandrel 
to  produce  yielding  (loading),  then  removing  the  mandrel  (unloading). 

The  region  where  yielding  occurs  during  expansion  is  determined  by  a 
loading  elastic-plastic  boundary,  p,  and  the  region  where  yielding  occurs 
during  unloading  is  determined  by  an  unloading  elastic-plastic  boundary, 
p'.  Assumptions  of  plane  stress  and  radial  symmetry  are  employed  in 
modeling  the  process.  Deformation  theory  is  used  to  relate  plastic 
strains  to  total  stress  and  the  Von-Mises  yield  criterion  is  used  to 
define  the  constant  yield  surface  for  plastic  flow.  The  continuum  equa¬ 
tions  are  averaged  through  the  thickness  of  the  plate. 

For  the  loading  process,  the  equations  consistent  with  these 
assumptions  are  (Reference  13). 
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where  the  variables  u,  n  ,  a„,  and  x  depend  on  r,  the  distance  from  the 

r  b 

center  of  the  hole,  and  represent,  respectively,  the  averaged  radial  dis¬ 
placement  after  the  loading,  the  radial  and  circumferential  components  of 
stress  after  loading,  and  the  proportionality  factor  associated  with  plas¬ 
tic  flow  during  loading.  Note  that  du/dr  and  u/r  are  the  radial  and 
circumferential  strains  —  the  first  terms  on  the  right-hand  side  of 
Equations  1  and  2  represent  elastic  strains  and  the  second  terms,  plastic 
strains.  The  material  constants  E'  and  v'  are  respectively  the  modified 
Young's  modulus  and  modified  Poisson's  ratio  (Reference  14)  given  by 


E' 


E(l+2v) 
2(1  +v ) 


and 


V 

1  +v 


where  E  and  v  are  the  usual  Young's  modulus  and  Poisson's  ratio.  The 
material  constant  aQ  is  the  yield  stress  determined  from  the  uniaxial 
stress/strain  curve  and  is  equal  to  the  effective  stress  for  perfectly 
plastic  modeling  (Reference  15). 

For  the  unloading  process,  the  constitutive  relations  (Equations  1 
and  2)  are  modified  to  include  the  permanent  deformations  due  to  loading 
according  to  Reference  16. 


37 -!■  (tr  ••'V*!1,  <w  *  Jx<V«e>  (5) 

w/r  -  1,  (te-v'tr)  +  1  »'(2te-tr)  ♦  1  »(2v«r)  (6) 

where  w,  tr>  tQ,  and  x',  represent  respectively,  the  averaged  residual 
radial  displacement  field,  averaged  radial  and  circumferential  components 
of  residual  stress,  and  proportionality  constant  associated  with  plastic 
flow  during  unloading.  Observe  that  the  residual  strains  dw  and  p- are 

comprised  of  elastic  strains  due  to  locked  in  residual  stresses  and 
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permanent  deformations  resulting  from  plastic  flow  during  loading  and 

unloading.  The  residual  stresses  and  unloading  plastic  proportional ity 

factor  satisfy  Equations  3  and  4  with  a.  o. ,  and  X  replaced  by  t  ,  t„, 

r  0  r  0 

and  X ' . 

It  is  not  transparent,  but  it  has  been  shown  in  References  13,  16 
that  the  relations  in  Equations  1-6  uniquely  determine  u,  a  ,  aQ,  X,  w, 
t  ,  t„,  and  X'  when  four  consistent  boundary  conditions  are  given.  For 
the  loading  process,  the  physical  boundary  conditions  are  the  mandrel 
interference,  uQ,  and  the  free  surface  at  the  outer  radius,  i.e., 

u(a)  =  uQ  op  (b)  =  0 

characterizes  loading.  Both  inner  and  outer  surfaces  are  free  in  the 
unloading  process;  thus 

tr(a)  =  0  tp(b)  =  0 

characterizes  unloading. 

It  should  be  emphasized  that  because  of  radial  symmetry,  deformation 
theory  as  applied  above  leads  to  the  same  results  as  the  Prandtl -Reuss 
flow  theory  (Reference  15).  The  Prandtl-Reuss  theory  is  the  basis  for 
the  finite  element  analysis  implemented  by  various  researchers  (References 
8,  10,  11,  12)  investigating  plane  stress  problems  similar  to  the  cold¬ 
working  problem.  Thus,  the  analysis  herein  is  a  semi-exact  solution  to 
compare  with  finite  element  analyses  except  that  some  of  the  finite  ele¬ 
ments  analyses  include  the  effect  of  strain  hardening.  A  strain  harden¬ 
ing  model,  e.g.,  a  Ramberg-Osgood  representation  of  the  uniaxial  stress/ 
strain  curve,  could  be  added  to  the  relations  (Equations  1-4)  at  a 
moderate  expense  in  the  computational  effort  required  to  solve  the 
equations . 
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2.  METHOD  OF  SOLUTION 

Given  geometry,  b/a,  and  material  constant,  E,  v,  o0>  a  solution  to 
the  nonstandard  two  point  boundary  value  problems  which  model  the  cold¬ 
working  process  is  uniquely  determined  by  specifying  uQ.  However,  it  is 
shown  in  Reference  13  that  uQ  must  satisfy  0  <  uq  <  u*,  where  u*  is 

a  maximum  radial  displacement  for  which  the  equations  have  a  solution. 

The  parameter  u*  is  interesting  as  a  candidate  for  an  optimum  interference 

level  and  is  discussed  later.  The  field  quantities  in  the  loading  process 

are  implicitly  determined  by  the  loading  elastic-plastic  boundary  --  the 

implicit  algebraic  relations  are  presented  in  the.  appendix.  Hence,  there 

is  an  algebraic  relation  between  uQ  and  p.  It  is  shown  in  Reference  13 

that  this  relation  is  one-to-one  for  u  which  satisfy  ue  <  u  <  u*  where 

o  o—o—o 

is  that  mandrel  interference  that  just  produces  yielding. 

The  procedure  used  to  solve  the  two  point  boundary  value  problem 
was  to  determinp  p  from  uQ  using  a  Newton-Raphson  scheme.  With  p  deter¬ 
mined,  o  ,  o. ,  and  u  are  expressed  analytical ly .  It  turns  out  that  the 
r  0 

unloading  elastic-plastic  boundary  p'  can  be  expressed  analytically  in 

terms  of  p.  Then  the  residual  fields  t  ,  tQ,  and  w  can  be  obtained  from 

r  b 

implicit  analytical  expressions  in  p  and  p'  except  that  the  unloading 
plastic  flow  proportionality  factor,  x‘,  must  be  determined  by  numerical 
integration.  Thus  that  portion  of  residual  displacement  field  where 
plastic  flow  occurs  during  unloading  is  the  only  field  quantity  that  is 
computed  numerically,  and  this  numerical  computation  is  quite  amenable  to 
accuracy  checks.  Computations  for  a  complete  solution  require  less  than 
a  second  of  CDC  6600  central  processor  time;  thus  parameter  studies  are 
quite  inexpensive. 

3.  PARAMETRIC  SOLUTIONS 

Mathematically,  there  are  four  situations  that  can  occur  depending 
on  the  size  of  the  manarel  interference.  They  are: 
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a.  purely  elastic  loading  and  unloading 

b.  plastic  flow  during  loading  only 

c.  plastic  flow  during  loading  and  unloading 

d.  no  solution 


The  interesting  situations  for  modeling  the  coldworking  process  are  b 
and  c  which  occur  when  uQ  is  in  the  interval  (uQ,  uQ).  In  order  to  con¬ 
cisely  summarize  results,  the  following  normalizations  are  introduced  for 
stress  and  displacement  fields 


a  = 


a 


a  ,!Cj3 

a% 


and  the  distances  r,  p,  and  p'  are  given  in  terms  of  hole  radii.  It  turns 
out  that  all  dimensionless  curves  depend  on  b/a  and  v  but  are  universal 
for  material  properties  E  and  aQ. 

General  results  for  the  elastic-plastic  equations  are  presented  in 
Figures  2-8.  Poisson's  ratio  was  fixed  at  v  =  0.3.  The  uq,  p  relation¬ 
ship  is  presented  in  Figure  2  for  values  of  b/a  ranging  from  b/a  =  3  to 
b/a  =  40.  As  discussed  in  the  appendix,  the  maximum  value  of  p  varies 
from  p  =  2.96  for  b/a  =  2.96  to  p  =  1.75  for  infinite  b/a.  These  curves 
are  universal  for  the  elastic-plastic  problem  except  for  a  small  depen¬ 
dence  on  v  which  will  be  discussed  later. 

Stress  and  displacement  fields  for  mandrel  in  and  mandrel  out  situa¬ 
tions  are  presented  in  Figures  3-6.  Dimensionless  mandrel  interference 
is  varied  over  the  range  (u®,  u*)  and  b/a  is  fixed  at  40  (essentially  an 
infinite  outer  radius).  In  Figure  3,  the  corners  in  the  hoop  stress  curves 
occur  at  the  loading  elastic-plastic  boundary,  p.  The  additional  corner 
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in  the  residual  hoop  stress  curves  occur  at  p' (Figure  5).  The  effects 
of  finite  geometry  on  stress  and  displacement  fields  are  illustrated  in 
Figure  7  --  the  parameter  b/a  is  varied  over  the  same  range  as  in  Figure 
2.  The  dimensionless  interference  is  fixed  at  7  corresponding  to  a 
situation  where,  according  to  Figure  2,  geometry  affects  p.  The  relation¬ 
ship  between  mandrel  interference,  loading  elastic-plastic  boundary,  and 
residual  deformation  at  the  hole,  wQ,  is  presented  in  Figure  8.  These 
data  are  limited  to  a  single  geometry,  b/a  =  40.  An  Interpretation  of 
the  data  displayed  in  Figures  2-8  is  presented  in  the  next  section. 

Poisson's  ratio  has  little  influence  on  the  dimensionless  solutions 
for  even  moderate  coldworking  levels;  hence  no  curves  are  presented  to 
illustrate  its  effect.  Poisson's  ratio  enters  the  problem  through  the 
elastic  strains;  so  a  quantitative  feel  for  the  influence  of  Poisson's 
ratio  is  obtained  by  taking  the  ratio  of  elastic  to  total  strains  at  the 
hole. 
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V 


SECTION  III 
DISCUSSION  OF  RESULTS 

1.  SIGNIFICANCE  OF  PRESENT  RESULTS 

It  is  generally  agreed  that  the  increase  in  fatigue  life  observed 
for  coldworked  holes  is  attributable  to  the  residual  compressive  hoop 
stresses  at  the  hole.  The  residual  hoop  stresses  tend  to  counteract  the 
tensile  stress  concentrations  which  would  otherwise  occur  in  remote  load¬ 
ing  situations.  Figures  3  and  5  indicate  that  the  loading  elastic-plastic 
boundary  is  a  good  measure  of  the  size  of  the  region  where  there  are  com¬ 
pressive  residual  hoop  stresses.  It  is  striking  how  the  region  of  compres¬ 
sive  stresses  is  enhanced  by  unloading.  This  suggests  that  interference 
fit  fasteners  might  be  effective  even  when  loosely  fitting  if  they  were 
initially  torqued  sufficiently  to  produce  substantial  yielding  at  the 
hole.  At  any  rate,  p  is  a  possible  measure  of  the  effectiveness  of  a 
given  mandrel  interference  in  suppressing  growth  of  radial  cracks;  hence, 
the  curves  in  Figure  2  relating  p  and  uQ  may  be  quite  useful  to  the  design¬ 
er  in  visualizing  the  coldworking  effects  for  different  mandrel  interfer¬ 
ences. 

j 

In  particular,  there  is  a  maximum  value,  pm,  for  the  loading  elastic- 
plastic  boundary  which  corresponds  to  u*,  the  maximum  allowable  mandrel 

Interference  (Figure  2).  Since  u*  leads  to  the  maximum  region  of  residual 

compressive  stresses,  it's  reasonable  to  interpret  u*  as  a  mathematical 

optimum  interference  level.  From  a  tolerance  standpoint,  it's  interest- 
ing  to  note  that  the  maximum  p  is  almost  obtained  for  a  broad  range  of 
mandrel  interferences  preceding  u*.  Experimentally  determined  optimum 

interference  levels  are  compared  with  u*  in  a  later  section. 
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In  view  of  the  minor  influence  of  Poisson's  ratio  on  solutions,  the 
data  presented  in  Figures  2-8  is  comprehensive  for  this  modeling  of  the 
coldworking  process  in  the  sense  that  the  effects  of  material  properties 
and  ring  geometry  are  displayed.  The  major  effects  of  material  properties 
enter  algebraicly  through  the  normalizations.  In  particular,  the  dimen¬ 
sionless  parameter  characterizes  the  primary  effect  of  material 

°o 

properties,  e.g.,  for  a  fixed  geometry,  the  optimum  percent  interference 


u*/a  scales  with  °o 
0  F73 


F7J 


The  major  effect  of  b/a  is  in  determining  the  largest  ring  of  yielding 
which  can  be  produced  by  coldworking.  Hence,  the  optimum  radial  expan¬ 
sion,  u$,  is  affected  by  geometry  (Figure  2).  Figure  7  indicates  that, 
for  fixed  uQ,  the  effect  of  b/a  on  residual  stress  field  near  the  hole 
is  insignificant,  and  both  Figures  2  and  7  indicate  that  infinite  geo¬ 
metry  is  a  good  approximation  for  b/a  >  10. 

The  uQ  and  p  vs  wq  curves  in  Figure  8  are  presented  to  show  the 
possibility  of  determining  the  amount  of  coldworking  (as  indicated  by  p 
or  uQ)  by  simply  measuring  the  hole  diameter  before  and  after  the  cold¬ 
working  process.  The  actual  mandrel  interference  may  not  be  easy  to 
measure  if  the  mandrel  deforms  as  well  as  the  material  around  the  hole. 
At  least  one  manufacturer  (Reference  17)  has  specified  the  coldworking 
process  in  terms  of  initial  and  residual  hole  diameters. 


2.  OPTIMUM  INTERFERENCE  LEVELS 

In  order  to  maximize  the  fatigue  life  improvement  due  to  the  cold¬ 
working  process,  the  designer  needs  to  select  the  proper  amount  of 
radial  expansion  uQ.  In  current  practice  the  optimum  interference  level 
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is  usually  based  on  a  series  of  fatigue  tests  (References  3,4).  Since 
the  optimum  expansion  would  be  expected  to  vary  with  material  and  hole 
diameter,  this  empirical  approach  can  be  both  excessively  expensive  and 
time  consuming  for  use  in  design. 

As  suggested  earlier,  the  beneficial  effects  of  coldworking  may  be 
characterized  by  the  loading  elastic-plastic  boundary  p,  since  p  deter¬ 
mines  the  size  of  the  residual  stress  fields.  In  particular,  there  is 
a  maximum  p  and  corresponding  maximum  uQ.  This  maximum  mandrel  interfer¬ 
ence,  u*,  may  be  considered  an  optimum  level  since  it  leads  to  the  largest 
region  of  compressive  hoop  stresses. 

For  convenience,  a  graphical  representation  between  uoE  ^  and 

CTo 

hole  radius  a  is  given  in  Figure  9  for  the  large  plate  case  (b/a=10). 
Poisson's  ratio  was  specified  as  0.3,  but  as  discussed,  varying  v  for 
other  practical  values  has  a  relatively  small  effect  on  the  results. 

Thus,  from  Figure  9,  one  may  specify  the  material  properties  E'  and 
0Q,  select  a  hole  radius,  and  readily  determine  the  radial  expansion  uq 

required  to  cause  the  maximum  elastic-plastic  boundary  p. 


In  order  to  test  the  hypothesis  that  u£  is  the  optimum  expansion 
from  a  fatigue  life  standpoint,  experimental  data  reported  in  References 
3  and  19  is  included  in  Figure  9.  The  data  of  Reference  3  were  obtained 
from  a  series  of  constant  amplitude  fatigue  tests  on  coldworked  holes 
in  which  the  amount  of  mandrel  Interference  was  the  variable.  The  speci¬ 
mens  contained  initially  unflawed  holes  which  were  drilled,  reamed,  and 
coldworked  in  a  manner  consistent  with  common  manufacturing  practice. 

The  holes  in  the  2024-T851  aluminum  specimens  were  unfilled  and  loaded 
remotely,  while  the  7075-T6  aluminum  members  were  subjected  to  partial 
load  transfer  through  a  fastener  placed  in  the  hole  as  would  occur  in 
an  actual  structure.  As  the  interference  level  was  Increased,  the 
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constant  amplitude  fatigue  life  (cycles  to  failure)  tended  to  approach 
a  maximum  value.  However,  the  increase  in  fatigue  life  gradually  levels 
off  as  the  interference  is  increased,  thus,  it  was  sometimes  difficult 
to  precisely  locate  the  optimum  interference  within  the  experimental 
scatter.  The  minimum  expansion  required  for  this  maximum  fatigue  life 
was  selected  as  the  optimum  interference  for  that  hole  size.  The  data 
point  on  Figure  9  from  Reference  19  is  the  result  of  a  similar  set  of 
experiments  with  interference  fit  fasteners  in  which  constant  amplitude 
fatigue  life  was  also  measured  as  a  function  of  fastener  interference. 

As  seen  in  Figure  7  the  experimental  data  generally  agree  quite  well 
with  the  hypothesis  that  u*  represents  an  optimum  value  of  radial  expan¬ 
sion  for  fatigue  life.  In  addition,  the  gradual  approach  to  the  experi¬ 
mental  optimum  is  consistent  with  the  flattening  of  the  p  vs  uQ  curve. 
Figure  2,  and  suggests  that  there  may  be  a  desirable  large  tolerance  in 
specifying  the  mandrel  interference.  In  this  regard,  note  from  Figure  5 
that  there  is  little  difference  in  residual  hoop  stresses  for  6.9  and  8.3 
dimensionless  interferences  and  8.33S  is  near  optimal.  A  key  point  here 
is  that  Figure  2  gives  the  designer  a  feel  for  near  optimal  as  well  as 
optimal  Interferences.  In  fact,  physical  limitations  of  actual  manufac¬ 
turing  equipment  may  preclude  reaching  the  theoretical  maximum  expansion 
values  in  some  cases.  The  mandrels  used  to  coldwork  the  steel  specimens 
sometimes  failed,  for  example,  at  the  loads  required  to  pull  the  mandrel 
at  the  larger  interference  levels  (Reference  3).  Thus,  both  the  relations 
between  p  and  uQ  shown  In  Figure  2  and  the  optimum  theoretical  values 
shown  In  Figure  9  should  be  quite  useful  to  the  designer,  since  they  pro¬ 
vide  a  fatigue  life  criteria  for  selecting  an  interference  level  for  a 
new  material  and  hole  size. 
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SECTION  IV 
CONCLUSIONS 

The  deformation  theory  of  plasticity  has  been  used  to  determine  the 
residual  stress  and  displacement  fields  caused  by  pre-expanding  a  hole 
in  a  circular  disk  with  an  oversize  mandrel.  Since  fatigue  rated  mechani¬ 
cal  fastening  systems  commonly  employ  interference  fit  or  pre-expansion  to 
develop  a  beneficial  residual  stress  field  about  the  fastener,  this  prob¬ 
lem  has  considerable  engineering  interest.  The  results  obtained  here  are 
presented  in  a  general  form  which  should  assist  in  predicting  the  perfor¬ 
mance  of  many  specific  fastener  applications.  In  addition,  by  explicitly 
using  radial  symmetry,  the  approach  is  inherently  more  accurate  and  econo¬ 
mical  than  finite  element  techniques.  The  method  could  be  extended  to 
include  strain  hardening  effects. 


In  Figure  3,  the  hoop  and  radial  stress  distributions  caused  by 
various  Interference  fits  are  presented  in  a  dimensionless  form  which 
allows  determination  of  the  stress  distribution  by  simply  inputting  the 
appropriate  material  properties  and  interference  conditions.  These  "man¬ 
drel  in"  curves  model  the  behavior  for  interference  fit  type  fasteners. 

The  residual  stresses  occurring  after  the  hole  Is  allowed  to  radially 
unload  are  given  In  a  simillar  general  form  in  Figure  5.  These  "mandrel 
out"  curves  represent  the  coldworking  process  which  Is  receiving  increased 
interest  for  use  In  fatigue  critical  areas.  Comparing  Figures  3  and  5, 
one  notices  that  thfe  radial  unloading  allows  a  much  larger  compressive 
residual  hoop  stress  to  be  developed  next  to  the  edge  of  the  hole  than 
the  corresponding  "mandrel  in"  configuration. 

An  analytical  optimum  radial  expansion  u*  is  proposed  which  would 
allow  the  engineer  to  maximize  the  fatigue  Improvement  In  the  coldwork¬ 
ing  process.  The  analysis  developed  here  indicates  that  there  is  a 
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maximum  loading  elastic-plastic  boundary  pm  developed  by  the  radial  expan¬ 
sion  which  Is  dependent  only  on  geometric  factors.  The  proposed  optimum 
interference  level  u£  is  the  radial  expansion  which  causes  p  .  For 

n:  ®  ITidX 

u*E  v3 

convenience,  0  is  plotted  versus  hole  size  in  Figure  9.  As  shown, 

CTo 

experimental  data  taken  from  the  literature  for  three  materials  agree 
quite  well  with  the  hypothesis  that  u*  Is  the  optimum  radial  expansion 
from  a  fatigue  life  standpoint.  Thus,  the  curve  of  Figure  9  should  allow 
the  designer  to  maximize  the  fatigue  Improvement  resulting  from  the  cold- 
working  process. 


14 


AFML-TR-79-4048 


REFERENCES 


1.  R.  J.  Gran,  et  al.,  Investigation  and  Analysis  Development  of  Early 
Life  Aircraft  Structural  Failures,  Technical  Report  AFFDL-TR-70-149 
Wright-Patterson  Air  Force  Base,  Ohio,  March  1971. 


2.  J.  A.  Regalbuto  and  0.  E.  Wheeler,  "Stress  Distribution  from  Inter¬ 
ference  Fits  and  Uniaxial  Tension,"  Experimental  Mechanics,  Vol .  10, 
No.  7,  pp  274-280,  July  1970. 


3.  J.  L.  Phillips,  Sleeve  Coldworking  Fastener  Holes,  Technical  Report 
AFML-TR-74-10,  Wright-Patterson  Air  Force  Base,  Ohio,  February  1974. 


4.  G.  J.  Petra k  and  R.  P.  Stewart,  "Retardation  of  Cracks  Emanating 
from  Fastener  Holes,"  Engineering  Fracture  Mechanics,  Vol.  6,  No. 2, 
September  1974,  pp  275-282. 


5.  A.  Nadai,  Theory  of  Flow  and  Fracture  of  Solids,  McGraw-Hill,  New 
York,  1950,  Chapter  33. 


6.  H.  T.  Jessop,  C.  Snell,  and  G.  S.  Holister,  "Photoelastic  Investigation 
on  Plates  with  Single  Interference-Fit  Pins  with  Load  Applied  to 
Plate  Only,"  The  Aeronautical  Quarterly,  Vol.  7,  November  1956, 
pp  297-314. 


7.  H.  T.  Jessop,  C.  Snell,  and  G.  S.  Holister,  "  Photoelastic  Investi¬ 
gation  on  Plates  with  Single  Interference-Fit  Pins  with  Load  Applied 
(a)  to  Pin  Only,  and  (b)  to  Pin  and  Plate  Simultaneously,"  The 
Aeronautical  Quarterly,  Vol.  9,  May  1958,  pp  147-163. 


8.  M.  Allen  and  J.  A.  Ellis,  Stress  and  Strain  Distribution  in  the 
Vicinity  of  Interference  FiT  Fasteners,  Technical  Report  AFFDL-TR- 
72-153,  Wright-Patterson  Air  Force  Base,  Ohio,  January  1973. 


9.  J.  H.  Crews  Jr.,  "An  Elastic  Analysis  of  Stresses  in  a  Uniaxially 
Loaded  Sheet  Containing  an  Interference-Fit  Bolt,"  NASA  Technical 
Note  NASA-TN-D-6955,  Langley  Research  Center,  Hampton,  Virginia, 
October  1972. 


10.  J.  H.  Crews  Jr.,  "An  Elastoplastic  Analysis  of  a  Uniaxially  Loaded 
Sheet  with  an  Interference-Fit  Bolt,"  NASA  Technical  Note  NASA-TN-D 
7748,  Langley  Research  Center,  Hampton,  Virginia,  October  1974. 


15 


AFML-TR- 7  9-4048 


REFERENCES  (CONTINUED) 


11.  W.  F.  Adler  and  D.  M.  Dupree,  Stress  Analysis  of  Coldworked  Fastener 
Holes,  Technical  Report  AFML-TR-74-44,  Wrlght-Patterson  AFB,  Ohio, 
March  1974. 


12.  D.  F.  Green  and  J.  P.  Gallagher,,  "Allowable  Stresses  for  Cracked 

Coldworked  Holes,"  Technical  Memorandum  AFFDL-TM-74-1 65-FBE,  Wright- 
Patterson  AFB,  Ohio,  August  1974. 


13.  T.  W.  Ting  and  R.  M.  Potter,  "The  Problem  of  Forced  Fittings," 

Archives  of  Rational  Mechanics  and  Analysis,  Vol .  58,  No.  1,  1975, 
pp  77-94. 


14.  A.  E.  H.  Love,  A  Treatise  on  the  Mathematical  Theory  of  Elasticity, 
4th  Ed.,  Dover,  New  York,  1944. 


15.  A.  Mendel  son.  Plasticity:  Theory  and  Application.  Macmillan  Co., 
New  York,  1968. 


16.  T.  W.  Ting  and  R.  M.  Potter,  "Problem  of  Forced  Fittings,  Part  II. 
The  Unloading  Process,"  Internal  Report,  Applied  Mathematics  Lab., 
ARL,  Wrlght-Patterson  AFB,  Ohio,  September  1973. 


17.  Technical  Literature  on  Seamless  Sleeve  Coldworking  Method.  J.  0  King, 
Inc.,  Atlanta,  Georgia  30318. 


18.  R.  M.  Richard  and  J.  R.  Blakelock,  "Finite  Element  Analysis  of 
Inelastic  Structures,"  AIAA  Journal ,  Vol.  7.  No.  3,  March  1969. 


19.  J.  H.  Crews,  Jr.,  "Analytical  and  Experimental  Investigation  of 
Fatigue  in  a  Sheet  Specimen  with  an  Interference-Fit  Bolt,"  NASA 
Technical  Note  NASA- TN-D- 7926,  Langley  Research  Center,  Hampton, 
Virginia,  July  1975. 


16 


AFML-TR-79-4048 

APPENDIX 

1.  ANALYTICAL  EXPRESSIONS  FOR  FIELDS 


Results  in  References  5,  13,  16  can  be  used  to  obtain  analytical 
expressions  for  stress  and  displacement  fields  after  loading  and  residual 
stress  field  after  unloading.  Stress  fields  after  loading  are  given  by 


sin  («l(r)+  tt/6) 


a  <  r  <  p 


o  [l+3(b/p  )  j  [l+(b/p)]p 


where  d  is  related  to  r  by 


and  0  is  related  to  p  through 


sin  (#J  =  l/[l+3(b/p )j 


The  corresponding  radial  displacement  field  after  loading  is  given  by 

C  2 E '  C(l-2v')a  (r)  +  (2a  (p)-a  (p))  exp  [/3(d  -0(r))]  a 

u(r)  =< 

(  '  pci  [o0(r)-v'a  (r)]  p  <  r  <  b 
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The  residual  stress  field  is  given  by 


tr(r) 


2o 

-  sin(^(r  )-tt/6) 

/3 


a  <  r  <  p 


(p')-o  fp1) 


1  - (b/p 1 ) 


[1  - (b/r )£]+  or(r)  p  <_  r  <_ 


t9(r)  - 


2  a 


sin  ( ( r )  +tt / 6 ) 


a  <  r  <  p 


tr(p'  )-or(p' ) 
1  -  ( b/ p  * ) 2 


[1 +( b/r )  ]  +  oQ(r)  p1  <  r; 


where  ^  and  r  are  related  by 
2  -tt7/3 

-a^/3  exp  6  =  exp  (-^(r)/!  cos  (^ ( r ) ) 

2r2 

and  p'  is  determined  from  Equations  A . 1 ,  A. 2,  A. 3,  and 


2 

b--=r--  [sin  (ip(p '  ) )  -  sin  («S(p'))]  =  [cos  (ip(p  * ) )  -  cos 

p  ' 

The  residual  displacement  field  is  given  by 

^t0(r)-a9(r)'v'  (tr(r)-crr(i"'  )+A‘  (r)(2tQ(r)- 

Ci^e (r)“CTe (r)-v'  (tr(r)-ar(r)]+u(r)  p'  : 
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.  b 

(A. 3) 

(«l(p'))] 

r(r)}+u(r)  a  <_  r  <_  p' 
r  <  b 
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where  X'(r)  satisfies  the  ordinary  differential  equation 

t1  ■  -[(2A'*3/E'4  <VV>T- 
subject  to  x 1 (p 1 )  =  0. 

2.  EXTREME  VALUES  FOR  p 

max 

The  maximum  loading  elastic-plastic  boundary  occurs  when  ar(a)  = 

-2a 

o  (the  maximum  stress  consistent  with  the  Von-Mises  yield  criterion). 
“73 

The  relation  between  {i  and  r  for  this  situation  is 

2  (f-) 

3  exP;> —  =  exp  (-fi(r)/3  cos  (0(r))  (A. 4) 

2/ 

Case  1  (Infinite  geometry):  From  Equation  A. 2,  infinite  b  implies  «Sp  = 

0;  then  Equation  A. 4  yields 

(P[nax/a)  =  [exp  U/3/3)/2]1/2  %  1.75 

which  is  the  smallest  possible  value  of  Pmax/a- 

Case  2  (Completely  plastic  ring):  It  can  be  shown  that  the  largest 

/  ma  x 

obtains  when  the  ring  is  completely  plastic,  i.e.,  where  pm,„  =  b. 

max 

The  largest  ratio  of  b/a  that  sustains  a  completely  plastic  ring 
is  given  in  Reference  5  as 

<b/a>L  *«e(If)/Ai 
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Figure  1.  Geometry  of  Coldworking  Process 
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b/a  =  DIMENSIONLESS  GEOMETRY  PARAMETER 


Dimensionless  Loading  Elastic-Plastic  Boundary, 
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DIMENSIONLESS  DISPLACEMENT 


Figure  4.  Dimensionless  Radial  Displacement  Field  (w  vs  r/a) 
After  Loading 
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u  «  DIMENSIONLESS  MANDREL  INTERFERENCE 
0 


DIMENSIONLESS  DISTANCE  FROM  CENTER  OF  HOLE  (r/a) 

Figure  5.  Dimensionless  Radial  Stress  After  Unloading  -  Mandrel 


8.31  u  =  DIMENSIONLESS  MANDREL  INTERFERENCE 
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Figure  7a.  Effects  of  Geometry  Parameter  (b/a)  on  Stress  After 
Loading 
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DIMENSIONLESS  DISTANCE  FROM  CENTER  OF  HOLE  (r/a) 


Figure  7b.  Effects  of  Geometry  Parameter  (b/a)  on  Residual  Stress 
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DIMENSIONLESS  RADIAL  DISPLACEMENT 


DIMENSIONLESS  ELASTIC/PLASTIC  BOUNDARY 
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HOLE  RADIUS  a  (CENTIMETERS) 

0  1.0  2.0 


Figure  9.  Comparison  of  Theoretical  and  Experimental  Values  for 
Optimum  Interference  Level  as  a  Function  of  Hole  Size 
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